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Abstract 

The trajectories of test particles moving in the gravitational field of a non-spherically symmetric 
mass distribution become affected by the presence of multipole moments. In the case of hyperbolic 
trajectories, the quadrupole moment of an oblate mass induces a displacement of the trajectory 
towards the mass source, an effect that can be interpreted as an additional acceleration directed 
towards the source. Although this additional acceleration is not constant, we perform a general rel- 
ativistic analysis in order to evaluate the possibility of explaining Pioneer's anomalous acceleration 
by means of the observed Solar quadrupole moment, within the range of accuracy of the observed 
anomalous acceleration. We conclude that the Solar quadrupole moment generates an acceleration 
which is of the same order of magnitude of Pioneer's constant acceleration only at distances of a 
few astronomical units. 
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I. INTRODUCTION 



Pioneer 10/11 are spacecraft which were launched more than thirty years ago to explore 
the outer region of the Solar system. They were also the first missions to enter the edge 
of interstellar space. Due to their especial structure, they also represent an ideal system to 
carry out precision celestial mechanics experiments which permit acceleration estimations 
to the level of 10~®cm/s^. After encountering the outer planets Jupiter and Saturn, the two 
spacecraft followed hyperbolic orbits near the plane of the ecliptic to opposite sides of the 
Solar system. A careful analysis of orbital data from Pioneer 10/11 has been reported in 3, ^ 
which indicates the existence of a very weak, long-range acceleration, ap = (8.74 ± 1.33) x 
10~^cm/s^, directed towards the Sun. The most intriguing result about this acceleration is 
that it is constant, with a single measurement accuracy averaged over 5 days. Moreover, 
it has been measured at distances contained within the approximate range of 20 to 60 
astronomical units (AU) from the Sun. 

The most conservative hypothesis about the origin of this anomalous acceleration is that 
this kind of effects on Pioneer's tracking data must be due to some technological reasons 
related to the configuration and equipment of the spacecraft. However, the analysis reported 
m seems to have ruled out all possible non-gravitational forces generated inside and 

outside the spacecraft. Moreover, the authors have taken into account the accepted values of 
the errors in the determination of the planetary ephemeris. Earth's orientation, precession, 
and nutation. More adventurous explanations of the anomalous acceleration suggest the 
existence of new physical phenomena, especially inspired by the fact that ap is of the order 
of cH, where c is the velocity of light in empty space and H is the Hubble constant. For 
instance, in P, U, 0] it is suggested that the expansion of the Universe induces local tidal 
forces which generate the anomalous acceleration. In [6] a Kaluza-Klein model is proposed 
in which the fifth dimension becomes dynamical with an expansion rate that generates ap. 
A scale dependent cosmological term has been used in [3| which affects the motion of test 
particles. A different explanation is suggested in ^ by means of a time-dependent New- 
tonian gravitational constant that produces a long-range acceleration. More sophisticated 
explanations have been proposed in the literature (see, for instance m , for a recent review) 
which include dark matter and stringy inspired scalar fields |9|, 110|| . 

In this work, we investigate the possibility of explaining the anomalous acceleration by 
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means of the Solar quadrupole moment. Clearly, the additional acceleration due to the 
quadrupole is not constant, but decreases with distance. In this sense the, quadrupole 
moment could immediately be ruled out as an explanation for ap. However, our approach is 
different. Since ap has been observed only at distances less than 60 AU, with an accuracy of 
around 15%, our goal is to determine whether the acceleration generated by the quadrupole 
is of the same order of ap and can be fitted within its range of accuracy. This idea is 
based upon the fact that the quadrupole moment of the Sun has been determined in an 
exact manner only recently and has been given the value of J2 = —2.28 10^^ ± 15%, 
where J2 is a dimensionless parameter to be introduced below. Before this value was known, 
different methods generated completely different values of J2, sometimes with differences of 
several orders of magnitude. (In fact, such a great ambiguity in the value of J2 could easily 
be used to explain ap in the complete range of observation!) With this new value, we will 
show that the quadrupole moment of the Sun produces an acceleration which is of the same 
order of ap at distances of a few All's, but is several orders of magnitude less than ap at 
distances between 20 and 60 All's. 

We will analyze the relativistic hyperbolic motion of test particles in the gravitational 
field of a mass with quadrupole moment. For the sake of simplicity, we will assume that 
the quadrupole is due only to the asphericities of the mass distribution, preserving the 
axial symmetry, and will neglect the rotation of the source. This allows us to use a static, 
axisymmetric solution of Einstein's vacuum field equations as the metric which describes 
the spacetime. In particular, we will use the Erez-Rosen metric for our analysis. Moreover, 
we will assume that the motion is confined to the equatorial plane of the source. This is 
in good agreement with the actual path of Pioneer 10/11 which move near the plane of the 
ecliptic. 



II. GRAVITATIONAL FIELD OF A MASS WITH QUADRUPOLE MOMENT 

In this section we present a static axisymmetric solution of Einstein's vacuum field 



12| . This solution generalizes 



equations which was first discovered by Erez and Rosen 
Schwarzschild solution to include an arbitrary mass quadrupole moment. Originally, this 
solution was obtained in prolate spheroidal coordinates, but for the purposes of this work 
spherical coordinates (t, r, 9, (p) are more appropriate. In these coordinates, it can be shown 
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that it is symmetric with respect to the equatorial plane, 9 = 7r/2, so that it becomes a 



geodesic plane. For the sake of simplicity we will limit ourselves to the investigation of 
motion on the equatorial plane. Accordingly, the Erez-Rosen solution can be written as 



ree 

ds^ = (^1 _ ^ j e^i^dt' - (^1 - ^) e^'i^^-'^Ur' - r\-^^^d^^ , (1) 

where the constant m represents the Schwarzschild mass, g is a constant parameter that 
determines the quadrupole moment, and ip and 7 are, in general, functions of r and 9. On 
the equatorial plane, these functions take the form 

1 
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The general form of the function 7 is rather cumbersome. Here we quote only the leading 
term in the limit r — 00 which is sufficient for the analysis we will perform in the next 
section. 

In general, the Erez-Rosen solution is asymptotically flat and this allows us to determine 
its multipole moments in a covariant manner. For instance, using the Geroch-Hansen [l^ 
definition of multipole moments one can show that the monopole is m and the quadrupole 
moment is (2/15)gm^. Higher multipoles are given in terms of the monopole and quadrupole 
moments. So the parameters m and q acquire a clear physical meaning. As expected, in the 
limiting case (? = 0, the Erez-Rosen metric (P) coincides with the Schwarzschild solution on 
the equatorial plane. The hypersurface r = 2m turns out to be singular, in accordance with 
the black hole uniqueness theorems. An additional singularity is situated at r — > 0. Outside 
the hypersurface r = 2m, the spacetime is completely regular. 

At large distances from the source, the expansion of the metric component gu leads to 
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Later on we will use this expression to determine the acceleration generated by the 
quadrupole moment. 



III. MOTION OF TEST PARTICLES 



Let us consider the motion of test particles on the equatorial plane of the Erez-Rosen 
metric (jT}. The geodesic equation can be reduced to a set of first order differential equations 
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due to the existence of two constants of motion pt and p^p associated with the time-translation 
and the axial symmetry, respectively. For the case of a massive test particle, the geodesic 
equations can be written as: 

-^)e^^^t = p,, (5) 
r2e-2«^V9 = p^ , (6) 



f2 = 



Vt 



* r / V f J f2 



,2 

(7) 



where a dot represents the derivative with respect to the affine parameter along the geodesic. 
It is possible to perform a detailed analysis of these set of geodesic equations jl^ and find 
out all the effects of the quadrupole moment on the motion of test particles. Here, however, 
we are interested only on hyperbolic motion. Hence, let us consider the special case r = r{ip) 
so that r = {dr /d(f)(f, and introduce the following notations: 

2m 2m , . 

u = , a = . (8) 

u p^ 

Then, the geodesic equation (|7j) reduces to 

{u'Y = F{u) := e-^"^ [v!" - + a'ue-^'^^ + a^g-^?^ (p^e'^'?'^ - l)] , (9) 

where a prime denotes the derivative with respect to the angle coordinate (p. Since the 
expression {u'Y is positive for a timelike geodesic, we must demand that the condition 
F{u) > be satisfied. Accordingly, from Eq.© we obtain 

aYt > (1 - u)e^'^'f' (a" + u\^'^^) . (10) 

Because the left-hand side of this inequality is a positive definite constant, this implies that 
the roots of the equation F{u) = must be contained in the interval (l—u) > or, according 
to Eq.(jHI), in the interval r G (2m, oo). As mentioned in the last section, the hypersurface 
r = 2m contains in general a curvature singularity. Thus, hyperbolic motion is confined to 
the region where the Erez-Rosen metric is regular. 

On the other hand, simple algebraic manipulations show that the equation F{u) = can 
be expressed as a quadratic algebraic equation for the function exp{2qip), i.e., 

e4.^ + ^e^^^ - 7^^^ = . (11) 

(1 — UjW' 
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Since exp{2qilj) > 0, we need to consider only the positive root of this equation which can 
be expressed as 

(12) 



2^ V^^^^T^"^ 



Now, the function ijj is positive in the interval r G (2m, oo). Then, from Eq.(|T^ we obtain 
that 



These conditions can be represented in an equivalent manner as 

f{u) := u^{u - 1) + a'^{u - 1) + a'^pf > or < for g > or g < , (15) 

respectively. On the other hand, as can be seen from Eq.(jHl), the function f{u) coincides with 
the function F{u) for vanishing quadrupole moment (g = 0), i.e., f{u) determines the motion 
of test particles in the Schwarzschild spacetime. This is a nontrivial result. In fact, suppose 
that the equation F{u) = has a root at m = mq- Then at uq the value of g is related to the 
value of Uq by means of (fT^ with u replaced by uq. If, for instance, g < 0, then according 
to Eq. p5|) the Schwarzschild function /(mq) is negative. In the case of hyperbolic motion, 
the value u = uq corresponds to a vanishing "velocity", i.e., u'{u = Uq) = or, equivalently, 
f = 0, a condition which can be satisfied only at the perihelion of the hyperbolic trajectory. 
At this point, the motion of test particles in the Schwarzschild spacetime is not allowed, 
because there f{uo) < 0. From the form of the functions F{u) and f{u) one can show that 
in general 

f{u) < F{u) for g < and f{u) > F{u) for g > . (16) 

Let us consider the case of an oblate configuration (g < 0). Hyperbolic motion corre- 
sponds to the case where F{u = 0) > and possesses (at least) one positive root, say U2. 
The motion is confined to the region < -u < ^2 with r2 = 2m/ U2 being the perihelion 
distance. At m = M2, according to (fTHIl . we have that f{u2) < and so no motion is al- 
lowed in the Schwarzschild spacetime. If we "move" from U2 towards -u = 0, we will find 
a point, say Ui, where f{ui) = 0. Starting from this point and up to the point u = 0, 
where f{u = 0) > 0, the motion is allowed in the Schwarzschild as well as in the Erez-Rosen 
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spacetime. For the trajectory in the Schwarzschild spacetime the perihehon distance is given 
by ri = Irajux and at this point, according to (fT^ . F{u = ui) > 0, that is, the trajectory 
in the Erez-Rosen spacetime would correspond to a non-zero radial velocity. So we see that 
the perihelion distance of a hyperbolic path becomes affected by the presence of an oblate 
quadrupole moment in such a way that 

r2 < ri , (17) 

a result that trivially follows from the fact that ui < U2- Thus an oblate quadrupole moment 
reduces the perihelion distance. 

Let us denote by rEni^p) ^^e hyperbolic path of a test particle in the Erez-Rosen spacetime 
with mass m and a negative quadrupole q, and by rs{if) the corresponding path in the 
Schwarzschild spacetime which is obtained from rEui^p) t>y setting g = 0. Notice that 
the mass m must be the same in order to be able to compare both trajectories. Clearly, 
rs{'^) must be a solution of the equation {u')^ = f{u), whereas ter^^) is a solution of 
{u'Y = F'lu) with the same mass. As we have shown above, at the perihelion we have that 
^ERifper) < i^sifper)- Morcovcr, siucc in this case f{u) < F{u) for all values of m, we obtain 
that in general 

TERi^) < rs{^) (18) 

at all points of the trajectory. The difference between these two trajectories will become 
smaller as r increases. Moreover, ter will approach rs asymptotically at r — > oo, where 
the contribution of the mass quadrupole moment becomes negligible. If we were analyzing 
the motion of a test particle in the gravitational field of a mass with negative quadrupole 
moment, ignoring the contribution of the quadrupole, we would expect a path (</?), whereas 
the measured path would be ter^'-p), which is always less than rs^fp). Then, one could 
interpret this result as due to an additional acceleration directed towards the source of the 
gravitational field. In fact, the presence of a quadrupole moment generates an "effective" 
acceleration which leads to deformations of the trajectories of test particles. In the next 
section we will estimate this effective acceleration. 

To finish this section we mention that in the case of a prolate quadrupole moment (g > 0) 
the effect is opposite, that is, the perihehon distance increases and the effective acceleration 
would be directed outwards the gravitational source. This case, however, is not expected to 
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occur in astrophysical bodies because any deviations from spherical symmetry are usually 
generated by rotation which leads to oblate configurations. 



IV. THE EFFECTIVE ACCELERATION 

To find the exact hyperbolic orbit of a test particle in the gravitational field described 
by the Erez-Rosen metric one must solve the differential equation © in the range where 
this kind of motion is allowed. This can be done by applying standard numerical methods. 
Another possibility is to consider the case of motion at large distances from the source, 
where we can use a Taylor expansion of the functions ip and 7 that enter the function F{u) 
in Eq.(jni), and preserve the leading terms only. The resulting function F{u) turns out to be 
a cubic polynomial of u. In this case, the differential equation Q can be solved analytically 
as V? = V^(w) (see for an analysis of the resulting solutions), an expression from which, in 
principle, we could derive the acceleration due to the quadrupole moment q. Nevertheless, 
the resulting expression is rather cumbersome and an analysis seems to be possible only by 
applying numerical procedures. 

To evaluate the acceleration induced by the quadrupole moment, we will use the 
analogy with the method applied to measure the Solar quadrupole. Indeed, for an axially 
symmetric distribution of matter, for instance the Sun, the contribution of the quadrupole 
moment to the exterior gravitational field can be expressed in the metric component gu as 



where J2 is a dimensionless constant, P2(cos^^) is the Legendre polynomial, and Rq is the ra- 
dius of the Sun. Evaluating this expression on the equatorial plane {9 = vr/2) and comparing 
it with Eq.(j3)), we can see that both expressions coincide if 



Then, from Eq. ()19j) we can read the expression for the gravitational potential whose deriva- 
tive with respect to the radial coordinate produces the leading term of the acceleration. 
Thus, we obtain 




(19) 




(20) 
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For the evaluation of this quantity we take, according to the value of J2 = —2.28 10~ 
for the Solar quadrupole moment and Rq ~ 6.96 lO^^cm for the photospheric radius of the 
Sun. Moreover, we express the distance r in astronomical units, i.e., r = n 1.496 lO^^cm, 
where is a positive real number. Thus, we obtain 

= -1 4.33 10~^cm/s^ . (22) 

So we see that at a distance of one astronomical unit, the order of magnitude of coincides 
with ap. Considering the uncertainty in the values of J2, -Ro and m, one can show that Og 
and ap are still of the same order of magnitude at distances smaller than n ^ 2.1. In the 
interval between 20 and 60 All's, where the anomalous acceleration has been measured, 
is several orders of magnitude less than ap. 

V. CONCLUSIONS 

In this work we have investigated the possibility of explaining Pioneer's anomalous accel- 
eration by means of the Solar quadrupole moment. In the analysis of the hyperbolic motion 
of a test particle in the gravitational field of a mass with a negative quadrupole moment, we 
have shown that the complete trajectory becomes closer to the gravitational source, an effect 
that can be interpreted as due to the presence of an additional acceleration, Og, generated 
by the quadrupole moment. 

Although Og is not a constant quantity, its value at distances very closed to one AU is of 
the same order of magnitude of ap. The discrepancy between Og and ap increases drastically 
with distance, indicating that it is not possible to explain the anomalous acceleration at 
distances where it has been measured for Pioneer 10/11 spacecraft. 
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